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In a recent Letter, Yoshimura and Arai [l[ claimed that the conventional phase stochastic differential equation 
O ! (SDE) used in @, S, || does not give a proper approximation to limit-cycle oscillators driven by noise, and proposed 
a modified phase SDE. Here we argue that their claim is not always correct; both SDEs are valid depending on the 
situation. 

£ — Since physical noise has an associated time scale and all oscillators have a characteristic rate of attraction, which 

of the two SDEs is appropriate depends on the relative sizes of these two scales. As a simple example, let us consider 
the Stuart-Landau (SL) model used in [l|, Q driven by a colored noise generated by the Ornstein-Uhlenbeck process 
(OUP) [|, which is rescaled such that the amplitude relaxation time explicitly appears while keeping the limit cycle 
and its isochrons invariant, 
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W{t) = {T -1 (l + ic) + iuj}W - T- l {\ + ic)\W\ 2 W + V2e£0), (1) 



where W is a complex variable representing the oscillator state, T is the relaxation time of the amplitude, c and u> 
are parameters, e is the noise intensity, and is OUP noise that is applied only to the real component of W for 
simplicity. is Gaussian-distributed, and its correlation function is given by (£(i)£(s)) = exp (— \t — s\/t) /(2t), 
which converges to 6(t) as r — ► 0. Thus, £(t) gives a colored-noise approximation to the Wiener process [5|. 

Introducing the amplitude R = \W\ and the isochron phase 4> = aigW — cln \W\ [3, Q, Eq. (P) can be written as 



R(t) = T- 1 (i?-i? 3 ) + v / 2ecos(0 + clni?) (2) 
(j>(t) =uj- V2eRT x {sin(0 + clnR) + ccos(0 + clnR)} £(t). (3) 



It is now clear that T actually determines the relaxation time of the amplitude R. The limit cycle in the absence of 
the noise (e = 0) is simply R(t) = 1 and <fi(t) = u>t + const. 

Two different SDEs have been previously derived describing this and other noisy oscillators. The non-agreement is 
due to the order in which the white- noise limit and the phase limit are taken [|| . The "conventional" model obtained 
; i ■ by taking the phase limit in the first has the form: 

a 



d<j){t) = [lo + eZ{4>)Z' {<!>)} dt + V2eZ((j))dw(t), (4) 

where the phase sensitivity (or response) function Z(4>) = — sin</> — ccos(f> in the present example. Yoshimura and 
Arai's modified phase model obtained by taking the white-noise limit first is given by 



dcj){t) = [ui + e {Z(4>)Z'(4>) + Y ((/>)}] dt + V2eZ((j))dw(t), (5) 

where the extra term Y(<f>) = (1 + c 2 ) sin(20)/2 

To see which of the two reduced phase SDEs ((H [5]) approximates the original noisy SL model Eqs. (J2j [3]) better, 
we compare the stationary phase probability density functions (PDFs) obtained by direct Langevin simulations of 
Eqs. (0 [3]) for different pairs of (r, T) with the PDFs obtained from the two phase SDEs ([U E|) by numerically solving 
the corresponding Fokkcr-Planck equations. We fix u> = 1, c = 2, e = 0.01, and vary r and T keeping tT — 0.001 
constant. 

Figure 1(a) shows the stationary phase PDFs obtained for two typical cases, T = 0.01 <C r = 0.1 and r = 
0.01 «T = 0.1. The conventional phase SDE ((4]) nicely fits the original model when T«t, whereas the modified 
phase SDE ^ is better when r <C T. Figure 1(b) shows mean-square errors of the approximate PDFs yielded by 
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FIG. 1: (a) Comparison of stationary phase PDFs obtained directly from Eqs. ([2] [3} with those obtained from the phase 
SDEs (HO for (T, r) = (0.1, 0.01) and (0.01, 0.1). (b) Mean-square errors of the approximate phase PDFs of SDEs (HE) from 
the original PDFs of Eqs. ([2] [3} plotted as functions of the amplitude relaxation time T(= 0.001/r). 



SDEs ((HE]) from the original PDF given by Eqs. (JHE1) as functions of T(= 0.001/r). It is clear that the conventional 
phase SDE Q gives a better approximation for T < T*, while the modified SDE (O is better for T > T*. 

Summarizing, we have demonstrated that the conventional phase SDE ([4]) is also a proper approximation to noisy 
limit cycles with sufficiently fast amplitude relaxation, which can be used as a starting point for further analysis. It 
is a natural generalization of the ordinary phase equation driven by smooth signals (which becomes evident when 
written in the Stratonovich Langevin form [5[), and it has a practical advantage of being completely determined by u) 
and Z ((/>), which are both experimentally measurable. The results in [l| are valid when literally white noise is given 
to limit-cycle oscillators from the outset. However, white noise is actually an idealization of physical processes with 
small but finite time correlations. One should be careful of competing small time scales involved in the problem when 
the white-noise limit is taken. For example, near the bifurcation of a limit cycle, the attraction to the limit cycle is 
slow so that the SDE ([5]) is reasonable. However, far from the bifurcation or for relaxation oscillators where the time 
scale T is very small, the conventional SDE (HJ) should be used. 

Finally, the above intuitive arguments can be made rigorous by using a multi-scale or projection-operator method [El, 

, which yields a family of effective phase SDEs depending on the ratio t/T [7|. It can be shown that the two 
SDEs (HI EJ) are actually two extreme cases corresponding to T — > 0, r > and T > 0, r — » 0, respectively. 

H. N. thanks Yoshiki Kuramoto and Kensuke Arai for useful comments. 



[1] K. Yoshimura and K. Arai, Phys. Rev. Lett. 101, 154101 (2008). 

[2] J. Teramae and D. Tanaka, Phys. Rev. Lett. 93, 204103 (2004); H. Nakao, K. Arai, and Y. Kawamura, Phys. Rev. Lett. 

98, 184101 (2007). 
[3] D. S. Goldobin and A. Pikovsky, Physica 351A, 126 (2005). 

[4] G. B. Ermentrout, R. F. Galan, and N. N. Urban, Phys. Rev. Lett. 99, 248103 (2007). 
[5] W. Horsthemke and R. Lefever, "Noise-Induced Transitions" (Springer, Berlin, 1984). 

[6] A similar situation arises in the analysis of classical Brownian particles with inertia driven by colored noises. See e.g. R. 

Kupferman, G. A. Pavliotis, and A. M. Stuart, Phys. Rev. E 70, 036120 (2004). 
[7] J. Teramae, H. Nakao, and G. B. Ermentrout, in preparation. 



